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Abstract. We compute deep-sea energy spectra and zenith- 
angle distributions of the atmospheric muons, both conven- 
tional and prompt. The prompt muon contribution to the 
muon flux underwater is calculated taking into considera- 
tion predictions of recent charm production models in which 
probed is the small- a; behavior of the gluon distribution in- 
side of a nucleon. We argue for a possibility to discriminate 
the PQCD models of the charm production differing in the 
slope of the gluon distribution, in measurements with neu- 
trino telescopes of the muon flux at energies 10-100 TeV. 



1 Introduction 

A correct treatment of the charm hadroproduction is impor- 
tant to the atmospheric muon and neutrino studies since short- 
Uved charmed particles, D^, D^, Df, A+, produced 
in collisions of cosmic rays with nuclei of the air, become 
the dominant source of atmospheric muons and neutrinos at 
energies ii^ ^ 100 TeV. Thus, one needs to take them into 
consideration as the background for extraterrestrial neutrinos 
(Learned and Mannheim, 2000). Muons originating from de- 
cay of these charmed hadrons are so called prompt muons 
(PM) that contribute to the total atmospheric muon flux. 

Another aspect of the interest to the charm production re- 
lates to the gluon density at small gluon momentum fraction 
X. The gluon density at small x is of considerable impor- 
tance because this influences strongly the charm production 
cross section, both total and inclusive. Recently Pasquali et 
al. (1999) and Gelmini et al. (2000, 2001) have analyzed the 
influence of small-a; behavior of the parton distribution func- 
tions (PDFs) to the atmospheric lepton fluxes at the sea level. 
Basing on next-to-leading order (NLO) calculations of the 
perturbative Quantum Chromodynamics (PQCD), they pre- 
dict PM fluxes at the sea level depending strongly on proton 
gluon distributions at small x scale, x < 10~^. 

The muon spectra underwater computed with the model 
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of Pasquali et al. (1999), in which used were the MRSD_ 
(Martin et al., 1993) and the CTEQ3M (Lai et al., 1995) 
sets of PDFs, were recently discussed (Naumov et al., 2000; 
Sinegovskaya and Sinegovsky, 2001). In this talk, using pre- 
dictions of more recent PQCD model (Gelmini et al., 2000, 
2001) of the charm production, we discuss the PM contribu- 
tion to muons underwater at depths typical for operating and 
constructing neutrino telescopes (AMANDA Collab., 2000; 
ANTARES CoUab., 2000; Baikal CoUab., 1997; NESTOR 
Collab., 2000). Namely, we try to study a PM flux under- 
water dependence on the slope A of the gluon PDF at small 
x: xg{x) oc The nature of the small-x behavior of 

the gluon density is now under extensive discussion (see, for 
example. Ball and Landshoff (2000); Schleper (2001); Vogt 
(2000); Yoshida (2001)). The small-x behavior of the PDFs 
is the subject of the deep interest since the reahzation of the 
underlying dynamics is yet far from being complete. 



2 PDFs and charm production models 

Due to dominant subprocess in heavy quarks hadroproduc- 
tion, gg — > cc, the charm production is sensitive to the gluon 
density at small x, where x is the gluon momentum fraction. 
One may evaluate the scale of x in cosmic ray interactions as 
follows. The product of the gluon momentum fraction xi of 
the projectile nucleon and that of the target x near the charm 
production threshold (~ 2mc) is xix = 4TO^/(2mAr£^o). 
where Eg is the primary nucleon energy in the lab frame. 
Since a muon takes away about 5% of the primary nucleon 
energy,_Eo — 20i?^, we have = Q.l{mc/mN){mc/ E^). 
Because of the steepness of the primary cosmic ray spectrum 
only large xi contribute sizeably to the atmospheric charm 
production, one needs adopt xi > 0.1. Taking ml ~ 2 
GeV^, one may find for i?^ > 100 TeV the range of impor- 
tance to be a; < 2 • 10^^. It should be stressed, this range 
is yet outside of the scope of the perturbative next-to-leading 
order global analysis of parton distributions (Martin et al., 
1999; Lai et al., 2000). 
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Fig. 1. Vertical sea-level muon flux predictions. (See Bugaev et al. 
(1998); Sinegovskaya and Sinegovsky (2001) for data references.) 

The spectral index A being used in the PQCD charm pro- 
duction models (Pasquali et al., 1999; Gelmini et al., 2001) 
covers wide range from about 0.1, the number that may be 
connected to the soft pomeron with intercept 1 + eo (where 
eo = 0.08), to about 0.5 in the 3-flavor scheme of the BFKL 
approach (Balitsky et al., 1978). 

In Fig. 1 shown are sea-level PM fluxes (added to the con- 
ventional muon flux) predicted with PQCD charm produc- 
tion models by Pasquali et al. (1999) (hereafter PRS) and by 
Gelmini et al. (2000, 2001) (GOV), which are used in our cal- 
culation of the deep-sea muon flux. Let us sketch out these 
models. 

PRS-1 : 

The PRS-1 model (dot lines in Figs. 1, 2) (identical with the 
PQCD-1 in Ref. (Sinegovskaya and Sinegovsky, 2001)) is 
based on the MRSD_ set (Martin et al., 1993). The PDF 
input parameters are the followings: xg{x, Qq) ^ as 
a; 0, 4-momentum transfer squared Qq = 4GeV^; the 
sea light quark asymmetry, u < c?, is taking into considera- 
tion^ the QCD scale in the minimal subtraction scheme (MS), 
A4'^ = 0.215 GeV, corresponds to the effective coupling at 
the Z boson mass scale as{M^) = 0.111. The factorization 
scale is hf — 2mc, the renormalization one is fiR = rric, 
where the charm quark mass, TOc, is chosen to be equal 1.3. 
The sea-level prompt muon flux was parameterized by au- 
thors (Pasquali et al., 1999) as 



= -5.9H 



s sr 

2 



0.290y + 0.143?;^ - 0.0147y 



(1) 



where y = \og^Q{E ^ / GeV) . 
PRS-2 : 

In the PRS-2 model (identical with the PQCD-2 in Ref. (Sine- 



govskaya and Sinegovsky, 2001)), CTEQ3M set (Lai et al., 
1995) was used. Corresponding inputs that were utilized in 
fliis model are Af^ = 0.239 GeV, as(A/|) = 0.112, 
rric — 1.3 GeV, ^xp — 2mc, hr = mc, and A = 0.286 
at Qq = 1.6 GeV. The corresponding PM flux was parame- 
terized as 



log,Q[Elcj,^''^^iE^)/icm-\-hT-'GeY^)] 



-5.79 + 0.345y + 0.105y^ - 0.0127y^ 



(2) 



PRS-3: 



In this model the CTEQ3M set was also used. Differing from 
PRS-2 in renormalization and factorization scales, iip = 
Hji = rric, this model appears the uncertainty due to the scale 
choice. The PM spectrum predicted is now 



log,Q[Ef,<l>^'^^{E,)/{cm-h-hT-^GeY')] 
0.156y2-0.0153y^ 



-5.37 + 0.0191y- 



(3) 



GGV: 



Here we present results for the model, among those discussed 
by Gelmini et al. (2001), which is based on MRST set of 
PDFs (Martin et al., 1999) with different values of the slope A 
intherangeO.1-0.5, > 1.25 GeV^; a^iMj) = 0.1175. 
The factorization and renormalization scales are: 



where 



TUT = (4 + "^c)^^^ = 1.25 GcV, 



and characteristic transverse momentum kx is of ~ rUc- 

In order to compute PM flux underwater we parameterize 
with Eq. (4) sea-level muon spectra of the GGV model (see 
Fig. 7 in Ref. (Gelmini et al., 2001)): 

0?'''^(SM)/(cm-2s-isr-iGeV-i) 

= A{E^IGeV)-^^°+^'y+^^y'+^'y'^ . (4) 

In Table 1 five sets of the parameters to Eq. (4) are presented 

Table 1. A dependent parameters of the prompt muon differential 
spectra at the sea level (Eq. (4)). 



A 


71,10-" 


70 


71 


72,10-" 


73, 10-' 


0.1 


3.12 


2.70 


-0.095 


1.49 


-0.2148 


0.2 


3.54 


2.71 


-0.082 


1.12 


-0.0285 


0.3 


1.80 


2.38 


0.045 


-0.82 


0.911 


0.4 


0.97 


2.09 


0.160 


-2.57 


1.749 


0.5 


0.58 


1.84 


0.257 


-4.05 


2.455 



for different values of the small- a; gluon PDF spectral index 
A. 

We calculated the conventional muon flux basing on the 
nuclear cascade model by Vail et al. (1986) (see also Naumov 
et al. (1998); Bugaev et al. (1998)). High-energy part of this 
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spectrum for the vertical may be approximated with formula 

(incm-^s-isr-iGcV-i) 



14.35^-3.672 Ei<E.^E2, 



at different zenith angles, we used PQCD PM fluxes calcu- 
lated only for the vertical direction, supposing the isotropic 
approximation for prompt muons to be a reliable at least for 
(5) 10"* < E^ < 10*^ GeV at zenith angles 9 < 80°. 



where Ei = 1.5878 x 10^ GeV, E2 = 4.1625 x 10^ GeV. 

Five thin lines in Fig. 1 present the sum of the conven- 
tional muon flux Eq. (5) and the GGV PM flux (Eq. (4)) 
corresponding to A = 0.1,0.2,0.3,0.4,0.5 (numbers near 
lines). Dot lines show the same for PRS models, Eqs. (1- 
3). For comparison there are also shown contributions due to 
the quark-gluon string model (QGSM) and the recombina- 
tion quark-parton one (RQPM) (Bugaev et al., 1989, 1998) 
(dash-dot and dash, respectively). Ratios of PM fluxes ac- 
cording above models to the conventional flux are shown in 
Fig. 2. As one can see, the cross energy of the PM fluxes and 
conventional one covers the wide region from 150 TeV to 
^ 3 Pe V, that is more than one order of the magnitude. 



3 Prompt muon component of the flux underwater 

Muon energy spectra and angle distributions of the flux un- 
derwater was computed with the method by Naumov et al. 
(1994). The collision integral in the kinetic equation in- 
cludes the energy loss of muons due to bremsstrahlung, direct 
e~^e~ pair production and photonuclear interactions. The 
ionization energy loss and the small-w part of the loss due 
to e+e~ pair production (t; < 2 • 10^'', where v is the frac- 
tion of the energy lost by the muon) were treated as contin- 
uous ones. In our calculations of underwater muon fluxes 



The prompt muon fraction of the flux underwater, Rpm, 
defined as ratio of the integral prompt muon spectrum to the 
conventional one, is presented in Fig. 3 for the depth of 4 km 
of the water equivalent (w. e.) and for cos6' — 0.2. As is 
seen from the figure, Rp„i related to the gluon density slope 
A = 0.5 is a factor 3 greater than that for A = 0.1 at > 10 
TeV. 

Zenith-angle distributions of the prompt muon contribu- 
tion at depths 1 ~ 4 km w. e., calculated for E^^ > 100 TeV, 
are shown in Fig. 4. Only GGV model predictions were used 
here with A = 0.1 (dash) and A = 0.5 (solid). For the ver- 
tical one can see Rp„i rising from about 0.2 at the depth of 
Baikal NT (1.15 km) to about 0.5 at the NESTOR depth (- 4 
km). For the larger zenith angles, 9 ~ 75°, this contribution 
becomes apparently sizable at depths 3 — 4 km. Differences 
in the predictions owing to a change of A, from 0.1 to 0.5 
(see h — 2 and 3 km w. e.), are also clearly visible: the ratio 
Rpm{^ — 0.5)/i?p„i(A = 0.1) at h = 2km w. e. grows 
from about 1.5 to about 5 as cos 9 changes from 1 to 0.2. 

Here we supposed no differences between the PRS and 
GGV models apart from those related to the charm produc- 
tion cross sections. Actually one needs to compare the pri- 
mary spectrum and composition, nucleon and meson pro- 
duction cross sections and other details of the atmospheric 
nuclear cascade being used in above computations. These 
sources of uncertainties would be considered elsewhere. 
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Fig. 4. Ratio of the prompt muon flux underwater to the conven- 
tional one as a function of cos 9. 



4 Summary 

In order to test the small-x gluon distribution effect we have 
computed deep-sea prompt muon fluxes using predictions of 
charm production models based on NLO calculations of the 
PQCD (Pasquali et al., 1999; Gelmini et al., 2001). The pos- 
sibility to discriminate the PQCD models, differing in the 
slope of the gluon distribution, seems to be achievable in 
measurements of the underwater muon flux at energies 10- 
100 TeV. 

Being hardly appeared at the sea level for energies up to 
10^ GeV (Figs. 1, 2), a dependence on the spectral index A 
of the small-x gluon distribution becomes more distinct at 
depths 3 — 4 km w. e. (Figs. 3, 4). At the depth of 4 km and 
at the angle of ~ 78° one could observe the PM flux to be 
equal, for A = 0.5, to the conventional one even for muon 
energy ^ 10 TeV (the cross energy). While for A = 0.1 the 
cross energy is about 70 TeV. For the high energy threshold, 
Ef^ > 100 TeV, and at ft, < 3 km w. e., the ratio Rpm is 
nearly isotropic up to ^ 60°. The "cross zenith angle" at a 
given depth, Oc{h), depends apparently on A: 

cos |a=o.5— 0-3, COS 6*^ |a=o.i— 0.1, h — 3kmw.e. 
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